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The energy-level structure of a quantum system plays a fundamental role in determining its behavior and 
manifests itself in a discrete absorption and emission spectrum 1 1]. Conventionally, spectra are probed via fre- 
quency spectroscopy 1 1, 2] whereby the frequency vofa harmonic driving field is varied to fulfill the conditions 
/S.E = hv, where the driving field is resonant with the level separation AE (h is Planck's constant). Although 
this technique has been successfully employed in a variety of physical systems, including natural Ol 0] and 
artificial (3ll4ll5ll^l7l[8ll9l[Tol[TTl ll2Lri3ll atoms and molecules, its application is not universally straightfor- 
ward, and becomes extremely challenging for frequencies in the range of lO's and lOO's of gigahertz. Here 
we demonstrate an alternative approach, whereby a harmonic driving field sweeps the atom through its energy- 
level avoided crossings at a fixed frequency, surmounting many of the limitations of the conventional approach. 
Spectroscopic information is obtained from the amplitude dependence of the system response. The resulting 
"spectroscopy diamonds" contain interference patterns and population inversion that serve as a fingerprint of 
the atom's spectrum. By analyzing these features, we determine the energy spectrum of a manifold of states 
with energies from 0.01 to 120 GHzxh in a superconducting artificial atom, using a driving frequency near 0.1 
GHz. This approach provides a means to manipulate and characterize systems over a broad bandwidth, using 
only a single driving frequency that may be orders of magnitude smaller than the energy scales being probed. 



Spectroscopy has historically been used to obtain a wide 
range of information on atomic and nuclear properties 
Early on, the determination of spectral lines helped eluci- 
date the principles of quantum mechanics through studies of 
the hydrogen atom and facilitated the discovery of electron 
and nuclear spin. Since then, several spectroscopy techniques 
to determine absolute transition frequencies (or, equivalently, 
wavelengths) have been developed, involving the emission, 
absorption, or scattering (e.g. Raman) of radiation. The ad- 
vent of tunable, coherent radiation sources at microwave and 
optical frequencies led to the age of modern atomic spec- 
troscopy, where a primary approach is to identify absorption 
spectra as the source frequency is varied [i, 2J ■ 

The study of artificial atoms at cryogenic temperatures, 
whose spectra extend into the microwave and millimeter wave 
regimes (10-300 GHz), faces different challenges. Stable tun- 
able microwave sources in excess of 70 GHz exist, but are 
expensive, and generally require multipliers which are inef- 
ficient and intrinsically noisy 1 14]. Frequency dependent dis- 
persion and attenuation, tight tolerances to control impedance, 
and multi-mode or restricted-bandwidth performance of trans- 
mission lines and waveguides ifTill . limit the application of 
broadband frequency spectroscopy in these systems. De- 
spite these challenges, superposition states in superconduct- 
ing 0, H Etl and semiconducting artificial atoms |0] have 
been probed directly up to several lO's of GHz. A number 
of leading groups have developed alternative techniques, e.g., 
resonant- and photon-assisted tunneling [st], which can be 
used to access spectroscopic information in specific systems at 
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even higher frequencies, though each has its own advantages 
and limitations and may not be easily applicable to other sys- 
tems. 

"Amplitude spectroscopy," introduced in this manuscript, 
probes the energy level structure of a quantum system via its 
response to driving-field amplitude rather than frequency (Fig. 
lA). It is applicable to systems with energy-level avoided 
crossings that can be traversed using an external control pa- 
rameter. For a generic artificial atom, such longitudinal ex- 
cursions throughout the energy level diagram (Fig. IC) are 
realized by strongly driving the system with an external field 
at a fixed frequency, which may be several orders of magni- 
tude lower than the energy-level spacing. In this limit, the 
system evolves adiabatically, except in the vicinity of energy- 
level avoided-crossings where Landau-Zener-type quantum 
coherent transitions occur The quantum interference between 
repeated Landau-Zener transitions gives rise to interference 
fringes that encode information about the system's energy 
spectrum. By trading amplitude for frequency, the amplitude 
spectroscopy approach allows one to probe quantum systems 
with strong coupling to external fields, such as solid-state arti- 
ficial atoms, over extraordinarily wide bandwidths, bypassing 
many of the limitations of a frequency-based approach. 

We demonstrate amplitude spectroscopy with a supercon- 
ducting qubit, a solid-state artificial atom exhibiting discrete 
energy states |7] that can be strongly coupled to external ra- 
dio frequency (RF) fields while preserving quantum coher- 
ence 19J. Artificial atoms are natural systems for probing a 
wide range of quantum effects: coherent superpositions of 
macroscopic states 12111,11], Rabi oscillations jol [Tol [TH, [121 
Us, 15, 16, 17, 1^, incoherent Landau-Zener transitions lll9ll . 
Stiickelberg oscillations ll20l l2ll l22i l23ll . microwave cool- 
ing Ii24l I25I I26I1 . and cavity quantum electrodynamics ll27l 
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FIG. 1: Amplitude spectroscopy with long-pulse driving towards saturation. (A) Amplitude spectroscopy diamonds. The qubit is driven at 
a fixed frequency v = 0.160 GHz, while amplitude V is swept for each static flux detuning 5/dc. Color scale: net qubit population in state 
\q, L), where L (R) labels diabatic states of the left (right) well of the qubit double-well potential, and q = 0,1, 2, ... labels the longitudinal 
modes. The diamond edges signify the driving amplitude V for each value of 5fdc when an avoided level crossing is first reached (amplitudes 
Vi — Vb for (5/dc = (J/dc)- The main diamond regions, symmetric about <5/dc = 0, are labeled Dl to D5. Arrows indicate signatures of 
transverse mode coupling (see Fig.|4l(. Top axis: the |0, L) — |0, R) energy spacing AEo.o accessed by V from 5/dc = 0. (B) Schematic of the 
qubit surrounded by a SQUID magnetometer readout. Static and radio-frequency (RF) fields control the state of the qubit: a 3-/is cooling-pulse 
followed by an amplitude-spectroscopy pulse of duration At and amplitude V. The qubit state is determined with a synchronous readout pulse 
applied to the SQUID. (C) Schematic energy-level diagram illustrating the relation between the driving amplitude V and the avoided crossing 
positions for a particular static flux detuning 5/dc ~ i5/dc- The arrows represent the amplitudes Vi — V5 of the RF field at which the illustrated 
avoided crossings are reached, marking the onset of the diamond regions in (A). 



l28i l29tl have been demonstrated with these systems. Signif- 
icant progress has also been made toward their application to 
quantum information science ifsoKsil] . including state initial- 
ization m, tunable IH SH^and long-distance iH HI] 
coupling, quantum control f37', '38",'3^,'4^, quantum state to- 
mography |41], and measurement 142. 43. 441. 

Our qubit (Fig.[T}3) is a niobium superconducting loop in- 
terrupted by three Josephson junctions (see also App.lAt ll45l 
l46(l . The qubit potential has a two-dimensional double-well 
profile near flux-bias / w $o/2, parameterized by the de- 



tuning (5/ = / — $0/2, where $0 is the superconducting 
flux quantum (Fig. A5). The system's energy eigenstates are 
comprised of transverse (p — 0,1,2....) and longitudinal 
(q — 0, 1, 2, . . .) modes, with energies controlled by the flux 
detuning Sf. When the potential is tilted so that resonant inter- 
well tunneling is suppressed, the eigenstates closely approxi- 
mate the diabatic well states localized in the left (L) and right 
(R) wells, and are associated with loop currents of opposing 
circulation. In this limit, the energies of localized states in 
the left (right) well increase (decrease) approximately linearly 
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with flux detuning (Fig. \Vp). Whenever the diabatic states 
\p, q, L) and \p' , q',R) are degenerate, resonant interwell tun- 
nehng mixes them and opens avoided crossings Apqyqi. 

We drive the qubit longitudinally with a time-dependent 
flux 6f{t) = (5/dc — $if sinojt, with a harmonic term of am- 
plitude $rf , proportional to the source voltage V, that induces 
sinusoidal excursions with frequency h' = u!/2Tr through 
the energy bands about a static flux bias (5/dc- Because the 
transverse and longitudinal modes should be decoupled for a 
symmetric system, we assume initially that only the lowest 
transverse mode is populated and use the reduced notation: 
\p,q,L) ^ |g,L), \p',q',R) ^ |g',R), and Apq,p,q, ^ A,,,, 
(Fig. A5B). We do observe, however, signatures of weak exci- 
tations of transverse modes (see Fig.|4]and related discussion). 

Each experiment uses the pulse sequence shown in Fig.[T}3, 
which consists of a harmonic cooling pulse to initialize the 
qubit to its ground state |24], followed by the desired ampli- 
tude spectroscopy pulse. The qubit state is determined with a 
synchronous readout pulse applied to a superconducting quan- 
tum interference device (SQUID) magnetometer (App. IA3l l. 
Using this technique, we investigate both the long and short 
time behavior of our qubit, and infer the energy-level slopes 
niq, and the splittings A, and locations Sfq^q' of avoided 
crossings to construct the energy level diagram. 

Fig. [T]\ displays the amplitude spectroscopy of the qubit 
driven towards saturation. We choose a driving frequency i/ 
such that, throughout the driving cycle, hi/ is generally much 
smaller than the instantaneous energy-level spacing, yet the 
speed is large enough to make the evolution through avoided 
crossings non-adiabatic. In this regime, Landau-Zener tran- 
sitions drive the system into coherent superpositions of en- 
ergy eigenstates associated with different wells. Four primary 
"spectroscopy diamonds" with large population contrast, cen- 
tered about (5/dc ~ (Dl, D2, D3, and D4), are observed 
in the data; they are flanked by eight fainter diamonds. The 
diamond structures result from the interplay between static 
flux detuning and driving amplitude, which determine when 
the various avoided crossings are reached. Because the on- 
set of each diamond is associated with transitions at a par- 
ticular avoided crossing, the diamond boundaries signify the 
avoided crossing locations. We use the linear relation between 

V and $rf (see Fig. |2]\) to obtain the avoided crossing loca- 
tions Sfq^q' listed in Table 1. 

For a particular static flux detuning (5/dc ~ Sf^^ < 
(Figs. [T}\ and [Tf) and driving amplitude increasing from 

V — Q, population transfer from |0,L) to |0,R) first occurs 
when the Ao,o avoided crossing is reached at y = Vi (front 
side of diamond Dl, see Fig. [T]\). For Vi < V < Yi, in- 
terferometric Landau-Zener physics |2ll, |22l 13 at the 
Ao,o avoided crossing results in the observed fringes (detail, 
Fig.|2]\). KvV — Vi, the adjacent avoided crossing Ai o is 
reached, inducing transitions between levels |0,R) and |1,L) 
and marking the back of D 1 . 

For Vi < V < V3, the saturated population depends on 
the competition between transitions at Ao.o and Ai o, on fast 
intrawell relaxation, and to a lesser extent on much slower 
interwell relaxation processes. Because in our experiment 
^0,0 hv K, Ai 0, the result is dominated by the dynam- 



ics at the Ai.o crossing. Although transitions |0,L) |0,R) 
are induced at the Ao.o crossing, stronger transitions at Ai o 
convert a substantial fraction of that population to 1 1 , L) . This 
excited population quickly relaxes back to |0,L), thus sup- 
pressing the net population transfer. For the combinations of 
(5/dc and V where interference between successive passages 
through A 1.0 is instead destructive, signatures of transitions 
due to Aq are visible, albeit with reduced contrast (detail, 

Fig.EK). 

At even larger amplitudes, transitions to additional excited 
states become possible. When V > V^, the qubit can make 
transitions between |0,L) and |1,R), marking the front side 
of diamond D2. The backside of this diamond is marked by 
the amplitude V = V4 that reaches A2.0, allowing transitions 
between |0,R) and |2,L). This description can be extended 
straightforwardly to the remainder of the spectrum. In this 
device, we did not find explicit signatures of coherent multi- 
path traversal between the Sf < and Sf>0 regions of 
the energy-level diagram (e.g., via avoided crossings Ai 1 and 

A2,2, ...)• 



There are several remarkable features associated with the 
amplitude spectroscopy shown in Fig.[T}\. First, we are able to 
probe the qubit continuously over extraordinarily wide band- 
width using a single driving frequency of only 0.16 GHz. The 
highest diamond (D5) in Fig. [T}\ results from transitions to 
energy levels more than 100 GHzx/i above the ground state. 
Even at such high energy levels, our artificial atom retains 
its energy-level structure in the presence of the strong driv- 
ing used to access them. 

Second, diamond D2 exhibits strong population inversion 
due to competition between transitions at avoided crossings 
Ao.i and Ai combined with fast intrawell relaxation to 
|0,L) and |0,R) (Fig. |2K). The transition rates at Ao,i and 
A 1.0 exhibit strong oscillatory behavior due to Landau-Zener 
interference, constructive or destructive, depending on the val- 
ues of 6fdc and V. The competition between these rates leads 
to a checkerboard pattern symmetric about (5/dc ~ with 
alternating regions of strong population inversion and effi- 
cient cooling. Similar checkerboard patterns are present in 
the higher diamonds. The population inversion observed here 
is incoherent, and can serve as the pump for a single-atom 
laser iH. 

The energy-level separation Ai?g^/ = h{\mq\ + 
l"^<?' l)('^/dc — ^fq.q') between states \q,L) and |g',R) is pro- 
portional to the net flux detuning from the location Sfq^q' of 
the avoided crossing A^ q/, and to the sum of the magnitudes 
of the energy-level slopes niq and m^/ . Because the relative 
phase accumulated between the \q,L) and [g', R) components 
of the wave function over repeated Landau-Zener transitions 
is sensitive to AEq^qi, the slopes can be derived from the in- 
terference patterns which arise when varying Sfdc- The A^th 
node in the interference pattern occurs when a relative phase 
of 2ttN is picked up during the qubit's excursion beyond the 
avoided crossing. The detuning locations of the nodes (inset, 
Fig.|2}J) follow the power-law Sqq^N'^/'^ with a prefactor Sqqi 
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FIG. 2: Energy-level slopes and interference patterns. (A) Detail of diamonds Dl and D2 (Fig. lA) showing interference patterns due to 
single (Dl) and multiple (D2) avoided crossings (see text). Diamond D2 exhibits strong population inversion. Arrows mark the location of 
the avoided crossing positions. Energy-level slopes are extracted in (B) from the interference fringes (dashed white lines) at 43 mVrms (Dl) 
and 150 mVrms (D2). The flux-to- voltage conversion factor is determined by the front side of Dl (dotted white line). (B) Determination of 
the energy-level slopes for levels |0, L), |0, R), \1, L), and |1, R). Detuning location of the TVth interference-node (see inset) vs. N^^'^ at 
the voltages marked with dashed lines in (A) and their corresponding linear fits (red line). Inset: Vertical slice from diamond Dl (43 mVrms)- 
Interference nodes TV used in the main panel are indicated by dotted lines. (C) 2D Fourier transform of both diamonds in (A). The sinusoids 
with wavenumbers fco and ki are contributions from diamonds Dl and D2, respectively, and are related to the energy-level slopes (see text). 



related to the energy-level slopes (App. lBTt : 




where a = 3tt/2 \/2. The factor a is the frequency-dependent 
conversion factor between RF flux and source voltage; its 
value a = 0.082 m<i>o/inVr,ns is inferred from the slope of 
the front edge of the first diamond (Fig. |2]\). Fig. |2J shows 
the N^^^ power-law fits to the nodes of the vertical slices 
in diamonds Dl and D2 which are used to extract toq 
mi (Fig. |2]\, dashed vertical lines), where we take \mq\ = 
\mq'\ = niq for q = q' in our system. The slopes are ob- 



tained sequentially from the fitted values Sqqi in Eq.\T\ starting 
with 2mo = 2.88 GHz/m<i>o, followed by mo + mi = 2.534 
GHz/m$o; their values are summarized in Table 1. 

The relation between the slopes niQ and mi is most clearly 
exhibited by the 2D Fourier transform of the diamonds in 
Fig. |2jZ! (App. IB 2l l. The observed structure in the first two 
diamonds arises from the underlying "Bessel-function stair- 
cases" of multi-photon resonances associated with transitions 
between the lowest four energy levels, where the n-photon 
absorption rate depends on driving amplitude through the 
square of the A^th-order Bessel function |2(||2l|2l. Us- 
ing Fourier analysis, the apparently complicated mesh of 
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FIG. 3: Amplitude spectroscopy with short-pulse driving. (A) Qubit response to a short RF pulse of variable length At as a function of static 
flux detuning Sfdc, with V = 181 mVims ; i' = 0.045 GHz. Top axis: driving field period (regions A-E); the maximum pulse width corresponds 
to '^l.S oscillation periods. The scan is performed at an amplitude value in the left side of diamond D3 (Fig.Al 1), which reaches all crossings 
through Ao,2 and A2,o. (B) Detail of the interference pattern in the boxed region of (A). The black curve marks the pulse width At at which 
the sinusoidal flux-excursion first exceeds and, subsequently, returns through Ao,2 at flux detunings S fdc- (C) Temporal oscillations along the 
horizontal line in (B) at the specific static flux bias Sfdc = Sfdc, fitted by a Landau-Zener model with damping (red line, see text). (D) A'^th 
interference node versus A'^^^^ along the vertical line in (B) and best linear fit. Inset: interference pattern along vertical line in (B) and node 
locations. 



overlapping Bessel functions is transformed to a pair of si- 
nusoids with periodicity related to the energy level slopes, 



ky = ±agsm (ksf^^/g), where g — 2(|m, 



The sinusoid associated with q = q' = arises from the tran- 
sitions at Ao,o, while the second sinusoid with q = 0, q' = 1 
and q = 1, q' — is degenerate and arises from the transi- 
tions at Aq.i and Ai o- All four diamonds and their individual 
Fourier transforms are presented in Figs. A6-A10. 

Valuable additional information about the energy level 
spectrum and temporal coherence is gained by perform- 
ing amplitude spectroscopy over short time scales (Fig. |3]l. 
Temporal-response measurements are performed by initializ- 
ing the system to the ground state at detuning 6fdc and then 
applying an RF field pulse of a variable length At, with fixed 
frequency and amplitude. The phase of the sinusoid is care- 
fully adjusted to maintain the timing and directionality of the 
RF-flux excursion through the energy levels between pulses. 



Crossing 

q,q' 


Location 

<5/q,q'(m$0) 


Coupling 


Slope 
m^' (GHz/m<E>o) 


0,0 





0.013 ±0.001 


1.44 ±0.01 


0,1 


8.4 ±0.2 


0.090 ± 0.005 


1.09 ±0.03 


0,2 


17.0 ±0.2 


0.40 ±0.01 


0.75 ± 0.04 


0,3 


25.8 ±0.4 


2.2 ±0.1 


0.49 ± 0.08 



TABLE I: Energy spectrum parameters of a superconducting artifi- 
cial atom determined using amplitude spectroscopy 



When the pulse ends abruptly at time At, the state of the sys- 
tem is preserved as the flux detuning instantaneously returns 
to Sfdc (finite decay- time corrections are discussed below). 

This is exemplified in Fig. [3] where parameters are tuned to 
investigate the A2,o level crossing (Fig. All). For example. 
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FIG. 4: Identification of transverse states of the qubit's double-well potential. (A) Qubit response to a short RF pulse of variable length At as 
a function of static flux detuning <5/dc, with V = 179 mVrms; i' = 0.025 MHz. The scan is performed at an amplitude value in the right side of 
diamond D3 (Fig. A12), where the crossing Aoa.oo (but not Aoo.os) is reached. The signatures of two crossings with transverse states, Ai2,oo 
and A22,oo, are indicated. (B) Schematic energy-level diagram showing the locations of the transverse states. 



starting in the ground state at positive detuning (Sf^c > 0), the 
qubit is driven through Sf{t) < 0, diabatically crossing Aq^o 
and Ai o at the beginning of the first quarter-period with no 
significant population transfer (beginning of region A). Sig- 
nificant population transfer first occurs in region A when A2,o 
is reached. The onset of population transfer is followed by 
brief temporal oscillations, which we use to find the splitting 
A2,o (see below). The population becomes stationary after 
the qubit returns through A2.0 in the second quarter period 
(region B). 

In the third quarter-period (region C), the driving pulse 
takes the qubit through the excited state avoided crossing Ao,2 
located on the opposite side of the level diagram (positive flux 
bias. Fig. [Tp). This event is marked by a second sharp pop- 
ulation transfer. The population subsequently remains nearly 
constant (region D) until a third abrupt population transfer oc- 
curs during the first quarter of the second period (region E), 
which signals the qubit's return through A2 and the repe- 
tition of the driving cycle. The population transfer does not 
reach its furthest extent in flux during the first half-period (as it 
does for the subsequent half-periods) because our pulse shape 
has slightly lower amplitude for times smaller than 5 ns. 

The observed response over short time scales is not sym- 
metric about Sfdc = 0. When starting in the ground state at 
static bias Sfdc — Sf^^ < 0, the system is drawn deep into 
the ground state during the first half -period, without inducing 
any transitions (see Fig[T]Z!). It is only during the second half- 
period that Ao,2 is finally reached and significant population 
transfer occurs. The detailed time dependence of population 
in this interval is shown in Fig.[3}i. The curved line in Fig.[3}j 
marks the pulse width At at which the sinusoidal flux excur- 
sion first exceeds and, subsequently, returns through Ao,2 for 
different flux detunings 6fdc', the sinusoidal excursion about 
the specific static flux Sf^^ is correspondingly indicated in 
Fig-TP- 

The temporal oscillations displayed in Fig.[3jJ and Fig.[3lZ! 



can be understood qualitatively as a Larmor-type precession, 
or "ringing," that results after the qubit is swept through the 
avoided crossing. In a pseudo-spin- 1/2 picture where the 
qubit states are identified with up and down spin states rela- 
tive to a fictitious z-axis, the qubit precesses about a tipped 
effective magnetic field which steadily increases in magni- 
tude and rotates toward the z-axis. This picture is consis- 
tent with a temporal analysis of the standard Landau-Zener 
problem, in which a linear ramp with velocity v sweeps the 
qubit through the avoided crossing. In the perturbative (non- 
adiabatic) hmit (IB 3b . this model yields the transition proba- 
bility 



Pit) 



^0,2 



(2) 



The integral in Eq. |2] often arises in the context of optical 
diffraction, giving rise to Fresnel oscillations similar to the 
coherent oscillations observed in Fig.[3f! . 

Although Eq. |2] captures the essential features of the data 
in Fig. [3f , to obtain a quantitative fit we must account for 
the non-abrupt ending of the pulse. This transient slightly 
modifies the total precession phase accumulated before read- 
out, and adds a small Mach-Zehnder-type interference due 
to the finite ramp speed back through the avoided crossing 
Ao,2- We found remarkable agreement between the data and 
a simulation of the Bloch dynamics of the two level system 
near Ao,2 which included longitudinal sinusoidal driving up to 
time t = At followed by a rapid decay over approximately 2 
ns, and intrawell relaxation with a rate of 0.65 ns^^ (Fig.[3]Z!). 
The value of Ao,2 is extracted as a fitting parameter and, in 
this regime, is largely insensitive to the details of the pulse 
decay and intrawell relaxation. 

As in the case of the long-time data, the energy-level slopes 
can be extracted from the vertical fringes (Fig. [3j3) using the 
A^^/-^ power-law fitting (Fig.[3jD) and Eq.[T] We used Eq.[T]to 
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infer m2 and from the sums toq + ^2 = 2.189 GHz/m$o 
and niQ + — 1.929 GHz/m$o. The short-time ampHtude 
spectroscopy procedure was applied to obtain Ag g/ for dia- 
monds D2-D4 and slopes rriq for diamonds D3-D4, and they 
are presented in Table 1 (Ann was previously obtained using 
the rate-equation approach 1221 )■ 

So far we have focused only on the strongly coupled lon- 
gitudinal modes. However, the lack of perfect symmetry al- 
lows us to probe excited transverse modes within our driving 
scheme as well. The population transfer is relatively weak, in- 
dicating small deviations from an ideally symmetric double- 
well potential and longitudinal driving. Signatures of these 
states appear in diamonds D3 and D4 (e.g., arrows in Fig.[T]\ 
and Fig. A13). The temporal response to a short RF pulse 
taken on the front side of diamond D3 (Fig. A12) is shown in 
Fig SK- The front side of diamond D3 results from accessing 
Ao2,oo during the first half-period, where some population is 
transferred from the right to the left wells (Fig. |4}\), and we 
have used the full notation explicitly indicating both longi- 
tudinal and transverse modes. Two weak population trans- 
fers are identified during the second half-period between the 
known positions of the longitudinal avoided crossings A00.02 
and Aoo.o3- This result is in agreement with simulations of 
the qubit Hamiltonian ||46|] . which indicate that two transverse 
modes |1, 2, R) and |2, 2, R) exist in this region, as illustrated 
in Fig.|4jj. Although we can identify their locations, the values 
of Aoo,i2 and Aoo,22 are not conclusively determined from 
this measurement, because the fringe contrast of their tempo- 
ral oscillations are small compared with those of the adjacent 
longitudinal crossings Aoo,o2 and Aoo,o3- 

The amplitude spectroscopy technique demonstrated here 
is generally applicable to systems with traversable avoided 
crossings, including both artificial and natural atomic sys- 
tems. Due to the sensitivity of interference conditions and 
transition probabilities to system parameters, it is a useful 
tool to study and manipulate quantum systems. The technique 
is extensible to anharmonic excitation, e.g., one can utilize 
arbitrary-waveform excursions through the energy levels and 
targeted harmonic excitations to achieve desired transitions. 
This type of hybrid driving has been very recently demon- 
strated with Cs |51] and Rb |52] atoms about Feshbach res- 
onances, systems containing weakly-coupled levels that are 
otherwise challenging to address with a standard frequency- 
based approach. The amplitude spectroscopy technique is 
complementary to frequency spectroscopy, allowing the char- 
acterization of a quantum system, and, more generally, should 
open new pathways for quantum control Issll . 
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APPENDIX A: SUPPLEMENTARY INFORMATION 

1. Device Fabrication and Parameters 

The device was fabricated at MIT Lincoln Laboratory. It 
has a critical current density Jc « 160A/cm^, and the 
characteristic Josephson and charging energies are Ej w 
(27r;i)300 GHz and Ec « (27rft)0.65 GHz respectively. 
The ratio of the qubit JJ areas is a « 0.84, and A « 
(27r/i)10MHz. The qubit loop area is 16 x 16 ^m^, and 
its self inductance is Lq « 30 pH. The SQUID Josephson 
junctions each have critical current Ico ~ 2 /lA. The SQUID 
loop area is 20 x 20 /im-^, and its self inductance is Ls ~ 30 
pH. The SQUID junctions were shunted with 2 1-pF on-chip 
capacitors. The mutual coupling between the qubit and the 
SQUID is M « 25 pH. 



2. Potential Energy of the PC Qubit 

The potential energy of the PC Qubit is a 2D anisotropic 
periodic potential with double-well structures at each lattice 
site. After designing for negligible inter-lattice-site tunneling, 
the potential energy can be visualized as a single double-well, 
as seen in Fig. A5a. It is convenient to parameterize the po- 
tential energy U with phase variables ipm — {'Pi — ^2)/'^ and 
(fip — {(pi + Lp2)l% where (p\ and 1^2 are the phases across 
junctions 1 and 2, respectively, in Fig. IB. It is also convenient 
to plot U/Ej, where Ej — <I>oIc/27r and /c is the critical current 
of the larger junction. When the double-well potential is sym- 
metric and the qubit is driven symmetrically, one can reduce 
the 2D potential to a ID double- well along the (/3,„ direction, 
as seen in Fig. A5b. This is the longitudinal direction along 
which the qubit circulating current varies through the phase 

a b 





i 







FIG. 5: 2D and ID double- well potential energies (see text), a, 

Contour plot of 2D double-well potential energy for our qubit at 
5/0 — 0.46"l'o, far detuned from the symmetry point Sfdc = 0. 
When the transverse quantum modes can be ignored, the potential 
energy can be treated as a ID double- well along the dashed line pic- 
tured, b, ID double-well potential energy for Sfo — 0.495"l'o- The 
wells are slightly tipped to the left and the four lowest energy eigen- 
states are shown. 
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3. Qubit Readout 

The qubit states are read out using a DC-SQUID, a sen- 
sitive magnetometer that distinguishes the flux generated by 
the qubit persistent currents, /q. The readout is performed by 
driving the SQUID with a 20-ns "sample" current /g followed 
by a 20-^s "hold" current (Fig. IB). The SQUID wiU switch 
to its normal state voltage Vg if Is > /sw,l (-^s > ^5w,r), 
when the qubit is in state \L) (|i?)). By sweeping the SQUID 
current /g and qubit flux detuning Sfo, while monitoring the 
presence of a SQUID voltage over many trials, we generate a 
cumulative switching-current distribution function. From this 
distribution, we extract a best-estimator line in the space of /g 
and Sfo that allows us to characterize the population of state 
\L) for a given flux detuning. 



4. Implementation 

The experiments were performed in a dilution refrigerator 
at a temperature of approximately 20 mK. The device was 
magnetically shielded with 4 Cryoperm- 10 cylinders and a su- 
perconducting enclosure. All electrical leads were attenuated 
and/or filtered to minimize noise. The electrical temperature 
of the device in the absence of microwave cooling was approx- 
imately 40 mK. After applying the microwave cooling pulse 
(Fig. IB), the effective temperature of the qubit was less than 
3 mK. Microwave cooling enabled the data to be obtained at a 
repetition rate of 10 kHz, much faster than the intrinsic equi- 
libration rate due to interwell relaxation. For all experiments, 
the static flux detuning was swept in 6-/i$o steps, and the RF 
amplitude was scanned in 0.5-mV steps (at the source). The 
pulse width was scanned in steps of 0.005 ns to .1 ns, and 
each data point comprised an average of 500 to 30,000 trials, 
depending on the desired resolution. 



Ao,o- Using the definition of the energy level slopes given 
in the text, this gives 

AE{t) « 2\mo\{aV - SU,) - \mo\aV uH"" , (B2) 

where we have fit the driving 5f{t) near the maximum of 
cos uot to a parabola. 

By setting AE{t^) = 0, we find the initial and final cross- 

ti.2 = Tt*, with U = In the 



mg times 

parabolic approximation to the driving signal, the phase ac- 
cumulated between crossings is 

4) ^ 2tt j {2\mQ\{aV - 6fdc) -\mo\aV ojh'^)i'^l,) 



Stt 

2\mo\ — {aV -6Uc)U 



(B4) 



Using the quantization condition for interference, (/> — 
2ttN, and the definition of t^,, we find the values of static 
flux detuning {<5/^^'} where interference occurs with driving 
source voltage V: 



2t:N ^ 2|r7io| 



3 [aVy/'^LO ' 

Rearranging Eq. (IB5b and using u; ~ 2ttv we find 

2/3 



(B5) 



N 



2/3 



(B6) 



where the prefactor to N"^/^ is identified with the fit parameter 
sq.o in the text. 

Expression (|B6l l can be generalized to any avoided crossing 
Ag.q' by making the replacement 2 1 mo I |TOg| + Irriq' |, from 
which we arrive at Eq. (1) in the text. 



APPENDIX B: SUPPLEMENTARY DISCUSSION 



1. Slope Extraction from Landau-Zener Interference Patterns 



2. 2D Fourier Transform of Amplitude Spectroscopy 
Diamonds 



The interference between sequential Landau-Zener transi- 
tions at an isolated avoided crossing is sensitive to the relative 
phase 

= 27r / AE{t')dt' (Bl) 
Jti 

accumulated by the two components of the wave function be- 
tween the first and second traversals of the avoided crossing. 
Here AE{t') is the instantaneous diabatic energy level sepa- 
ration at time t', and ii 2 are the times of the first and second 
traversals, respectively. We note that the energy AE is mea- 
sured in frequency units (GHz), and this is why the expression 
contains the factor 27r rather than 1 /h. 

For demonstration, we focus on the interference in the first 
diamond where the driving 5f{t) = —Sfdc + aVcosujt 
sweeps the qubit through only the lowest avoided crossing 



The amplitude spectroscopy plots in Figs. 1 and 2 of the 
main text display structure on several scales. On the largest 
scale, the boundaries of the "spectroscopy diamonds" are 
readily identifiable. The interiors of the diamonds are tex- 
tured by fringes arising from the interference between succes- 
sive Landau-Zener transitions at a single or multiple avoided 
crossings. On an even smaller scale, these fringes are com- 
posed of a series of horizontal multiphoton resonance lines. In 
order to extract information from these small scale structures, 
it is helpful to apply a transformation that is able to invert 
length scales; the two-dimensional Fourier transform (2DFT) 
provides this service. 

An analytic treatment is presented in detail in Ref. i50il : 
the main conclusions are presented here. The treatment in 
Ref. llsoll is applicable to the perturbative driving regime 
A'^/hv <C 1, where A is the splitting at the largest traversed 
avoided crossing and v — dE/dt is the speed of sweeping the 
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qubit through level crossing. In our device for the driving fre- 
quencies employed, this condition is satisfied in the first two 
diamonds, where we find good agreement between the analyt- 
ical treatment, numerics, and the data. For higher diamonds 
where the dynamics are non-perturbative, more complicated 
behavior is observed. The numerical approach can still be 
employed, but in practice we find that in such cases it is more 
efficient to extract the desired information directly from the 
short-time dynamics as described in the discussion of Fig. 3 
in the text. 

The internal structure of the first diamond arises from re- 
peated passages through a single weak avoided crossing. As 
discussed in |50], the primary features of the 2D Fourier 
Transform of qubit magnetization for this perturbative driv- 
ing regime can be understood by studying the 2DFT of the 
transition rate at the qubit level crossing. The 2DFT of the 
transition rate displays intensity concentrated along the curve 



fcf> = ± — sin 



(2 ^«/dc) ' 



(B7) 



where the flux detuning and the driving signal are measured 
in the energy units: (5/dc = 2|mo|(5/dc and V — 2\rao\aV . 
After going back to the physical units, Eq.( IB7l ) gives 



4|mo|Q; . 
kv — ± sm 



4|mo| 



■k 



Sfdc 



(B8) 



The simplicity of the result can be traced to the fact that the 
curve ( IB7I) reproduces time evolution of the quantum phase of 
the qubit |50] which is harmonic for harmonic driving. 

Most strikingly, the apparently distinct phenomena of inter- 
ference fringes and multiphoton resonances observed in the 
real space image are manifested as a single smooth curve in 
Fourier space. This structure can be understood by consider- 
ing kv and ks /^^ to be smoothly varying functions of the spa- 
tial coordinates (V, Sfdc)- Through a stationary phase analysis 
of the Fourier integrals, one finds the mapping between real 
space patches and regions of fc-space depicted in Fig. AlOa. 

In numerical simulations we found that the steady-state 
magnetization in the second diamond was well reproduced by 
a simple rate model based on incoherently adding two addi- 
tional transition rates to account for transitions at the avoided 
crossings with the left and right first excited states. These ad- 
ditional rates are calculated using values of A appropriate for 
the excited state avoided crossings (approximately 90 MHz), 
and also take into account the different slopes of the ground 
and excited state energy levels. 

In the above model, the Fourier image of diamond two is the 
sum of the Fourier transforms of the three relevant transition 
rates. Due to the difference in the dispersion of the lowest 
and second lowest qubit energy levels versus dc flux bias, for 
the Aoi and Aio level crossings we have (5/dc = (l™o| + 
|"ii|)(5/dc, V — (|™o| + |r7ii|)aV^. This gives a sinusoid 



^ 2(|mo| + |mi|)Q . 
kv = ± sm 



2(|mo| + |mi|) 



ksfdc (B9) 



with the wavelength different from that of ( IB8I 1 by the ratio 
of the slopes 2|mo|/(|mo| + Both sinusoids can be 



seen in the Fourier transform of the second diamond, shown in 
Fig. A7, which indicates that the transitions at the Aqo cross- 
ing and the Aqi (Aqi) crossings contribute to the observed 
pattern. From the measured ratio of the sinusoids' wave- 
lengths we obtain the ratio of the slopes of the qubit energy 
levels mo /mi without any fitting parameters. 

In addition, in the Fourier transform of the second diamond 
the real-space integration samples a more limited sector of the 
fringes arising from each of the avoided crossings as com- 
pared to the case of the first diamond (see Fig. AlOb). As 
a result, the regions around the middle part of a sinusoid, 
2( | mi|+|mo|) ^'^/dc ~ (" + l/2)7r, are absent from the sinu- 
soids arising from the Aqi and Aiq crossings, while the re- 

riTT, are absent from the 



gions around the nodes. 



4 1 mo I ''/d. 



sinusoid arising from transitions at the Aqo crossing. 



3. Perturbative Treatment of Fresnel Oscillations 

The time-dependent oscillations observed in temporal- 
response measurements (see Fig. 3C) result from Larmor pre- 
cession about a tilted axis following the qubit's transit through 
an avoided crossing. In the regime where the Landau-Zener 
transition probability is small, we use a perturbative model to 
relate these oscillations to the well-known Fresnel integral. 

By linearizing the sinusoidal driving signal near the mo- 
ment of traversal through the avoided crossing, t,, we ar- 
rive at the familiar Landau-Zener Hamiltonian H{t) = 
(ft/2) (ut(7'' + A(T^), where u = ft(|TOq| + jm,' |)$rf coswi* 
is the sweep velocity and A is the splitting at the avoided 
crossing. Next, we transform to a non-uniformly rotat- 
ing frame by |V'fl(t)) = e^-^^*)*' | V(0 where (t){t) = 
— ^J^vt'dt' — —jvP. The rotating frame Hamiltonian is 
purely off diagonal, with Aij {t) = A exp(--TOt^/2), 



Hr = 



ft 



Anit) 



2 V A|j(t) 0. 



(BIO) 



We now expand the system's time evolution operator 
U{t, to) to first order in A: 

C/(t,to)-l-- / HR{t')dt' + 0{A^). (Bll) 

This approach is valid when the driving conditions are far 
from adiabaticity, i.e. /v <C 1. The probability P{t) — 
\{1 \U{t,to)\ i)P to find the system in the excited state 
I t ) at time t given that it started in the ground state | | ) 
at vto ^ — A is given by 



A2 

Pit) = - 



to 



(B12) 



as given in Eq. (2) of the main text. The oscillatory depen- 
dence of this this integral on the final time t can be verified 
with the help of the Cornu spiral. 
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FIG. 6: Diamond 1 and its 2D Fourier transform, a, Diamond 1 
in Fig.l in the main text. Driving frequency v — 160 MHz. b, 2D 
Fourier transform of diamond 1. A single sinusoid is visible, along 
with its lower harmonics i50A . 



FIG. 8: Diamond 3 and its 2D Fourier transform, a, Diamond 
3 in Fig.l in the main text. Driving frequency u = 160 MHz. b, 
2D Fourier transform of diamond 3. Because the Fourier transform 
samples progressively smaller sectors of the diamonds as diamond 
number increases, the extent of the sinusoids in the kg direction is 
hmited Isoll . 




FIG. 7: Diamond 2 and its 2D Fourier transform, a, Diamond 
2 in Fig.l in the main text. Driving frequency u = 160 MHz. b, 
2D Fourier transform of diamond 2. Two sinusoids are visible with 
slightly different period due to the multiple energy bands and avoided 
crossings that lead to diamond 2 I50l1 . 



FIG. 9: Diamond 4 and its 2D Fourier transform, a, Diamond 
4 in Fig.l in the main text. Driving frequency f = 160 MHz. b, 
2D Fourier transform of diamond 4. Because the Fourier transform 
samples progressively smaller sectors of the diamonds as diamond 
number increases, the extent of the sinusoids in the kg /^^ direction is 
hmited {5^. 
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RF Source Voltage V (mV^^) 





RF Source Voltage V (nVfiyig) 



FIG. 12; Diamond 3 at = 25 MHz. Dashed line indicates the 
amplitude (179mV,ms) at which the pulse width scan in Fig. 4 was 
taken. 



a 



FIG. 10: Graphical interpretation of 2D Fourier transform tecli- 
nique. a, Diamond I at u — 400 MHz is pictured alongside a 
schematic of its Fourier transform. On a local scale within the 
wedge-shaped region of (V^, (5/dc) space outlined by a solid red line, 
the image everywhere looks like a simple series of evenly spaced 
horizontal bands. The Fourier transform over this region maps to a 
region of (ky, kg /^^ ) space localized near the kg axis as indicated 
by the red dot in (B). Within the region outlined by the dashed white 
line, the local periodicity is along the angled interference fringes; 
the Fourier transforms maps this region to the region of {ky, kgf^^) 
space localized near the extrema of the sinusoid in the kv direction 
as indicated by the open white circle, b, The Fourier integral samples 
a smaller sector of the real space image associated with the higher ex- 
cited states in higher diamonds. Due to the mapping between sectors 
of real space and localized regions of fc-space, the sinusoids associ- 
ated with higher diamonds are not fully developed. Only the portion 
of the sinusoid corresponding to the mapped region in real space is 
realized in the Fourier image. 




160 200 
RF Source Voltage V (mVrms) 

FIG. 1 1 : Diamond 3 at = 45 MHz. Dashed line indicates the 
amplitude (ISlmVrms) at which the pulse- width scan in Fig. 3 was 
taken. 



O 4 



3 

2 "2 
M 




160 



200 



240 



RF Source Voltage V (mVrms) 




TJ 

■a 



10 20 30 

RF Pulse Width tn (ns) 

FIG. 13: Pulse width scan at = 25 MHz in diamond 3. a, 

Diamond 3 at = 25 MHz. Dashed line indicates the amplitude 
(183mVims) at which the pulse width scan in Fig.A13b was taken, b, 
Pulse width scan in diamond 3 at i/ = 25 MHz. A loss of population 
occurs where the crossings Aoo,o3, A22,oo and Ai2,oo are reached. 
This is in contrast to Fig.4A, where crossing Aoa.oo was reached, but 
never Aqo.os. Also noted are the avoided crossings with which the 
diagonal lines in the diamond (see a) are attributed. 
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